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ABSTRACT 


We  study  three  eigenvalue  problems  of  boundary  layer  theory  first 
introduced  by  Fox  and  Libby,  *  and  pursued  further  by  Kotorynski.  ^  By 
completing  Kotorynski's  approximate  analysis,  we  exhibit  approximate 
eigenvalues  and  norms.  With  these,  the  convergence  of  the  eigenfunction 
series  solutions  of  two  heat  transfer  problems  can  be  investigated:  (a)  Flat 
plate  with  one  constant  wall  temperature  region  followed  by  another,  (b)  Flat 
plate  with  a  constant  wall  temperature  region  followed  by  an  insulated  wall. 
The  series  for  problem  (b)  is  found  to  converge  everywhere,  including  at 
the  station  of  change  in  the  boundary  condition.  The  series  for  problem  (a) 
does  not  converge  at  that  station,  and  convergence  for  both  appears  slow 
near  that  station.  Then  we  improve  Kotorynski's  approximation  and  use 
the  WKB  method  to  obtain  simple  formulas  for  the  eigenvalues,  which  prove 
very  accurate,  and  enable  the  eigenfunctions  to  be  found  by  a  simple  forward 
integration.  Finally,  we  study  problem  (b)  in  more  detail,  and  obtain  an 
exact  solution  by  a  numerical  finite  difference  procedure.  This  shows  the 
eigenfunction  series  solution  to  be  very  slowly  convergent  ne ar  the  beginning 
of  the  insulated  region.  The  solution  worked  out  by  Durgin°’  based  on 
Lighthill's^  approximate  method,  is  very  accurate  there.  In  addition,  it  is 
within  10%  of  the  numerical  solution  everywhere.  It  is  undoubtedly  the  most 
convenient  representation  of  the  solution,  agreeing  both  with  the  new  numeri¬ 
cal  solution,  and  with  Durgin's  experiments. 
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I.  INTRODUCTION 


Fox  and  Libby*  have  developed  perturbation  solutions  to  the  momen¬ 
tum  and  energy  equations  of  laminar  boundary  layer  theory.  The  momentum 
solutions  are  perturbations  of  the  Blasius  flat  plate  solution.  Among  the 
energy  solutions  are  those  for  initial  value  problems  with  the  Blasius  veloc¬ 
ity  field  and  either  constant  wall  enthalpy,  or  zero  heat  transfer,  downstream 
of  the  initial  station. 

In  their  investigations,  Fox  and  Libby  formulated  three  eigenvalue 
problems:  (1)  For  the  perturbation  of  the  momentum  equation  caused  by 
initial  values  differing  from  the  Blasius  profile;  (2)  for  the  initial  value 
energy  problem  with  constant  wall  enthalpy;  (3)  for  the  initial  value  energy 
problem  with  zero  heat  transfer.  In  each  case  they  found  the  first  ten  eigen¬ 
values,  norms,  and  eigenfunctions  by  a  numerical  searching  procedure. 

2 

Subsequently  Libby  recalculated  the  eigenvalues  and  norms  for  problem  (1) 

and  added  the  next  ten  sets  of  values. 

3 

Kotorynski  has  considered  these  three  eigenvalue  problems  analyti¬ 
cally  as  irregular  Sturm -Liouville  systems.  He  reduced  them  to  normal 
form  and  approximated  the  coefficient  function  by  the  sum  of  its  principal 
terms  near  the  wall,  and  far  from  the  wall.  He  pointed  out  that  the  resulting 
equations  could  be  solved  in  terms  of  Laguerre  polynomials,  and  found  approx¬ 
imate  expressions  for  the  eigenvalues.  The  resulting  eigenvalues  are  not 
particularly  accurate,  the  twentieth  for  problem  (1)  being  4.  3%  higher  than 
Libby's  computed  value.  Kotorynski  did  not  use  the  boundary  conditions  to 


find  explicitly  the  approximate  eigenfunctions  or  the  norms,  to  compare 
with  Fox  and  Libby's  numerical  values. 

The  purpose  of  the  present  paper  is  three-fold:  (1)  To  extend  the 
approximate  eigenvalue  results  of  Kotorynski  to  include  the  approximate 
eigenfunctions  and  norms,  and  investigate  the  convergence  of  some  solu¬ 
tions  obtained  by  Fox  and  Libby,  using  these  eigenvalues  and  norms. 

(2)  To  apply  the  results  of  a  WKB  approximation  to  obtain  another  approxi¬ 
mate  formula  for  the  eigenvalues  which  is  extremely  accurate  indeed,  and 
should  enable  numerical  calculation  of  the  eigenfunctions  and  norms  with 
little  difficulty  in  a  minimum  of  computer  time.  (3)  To  discuss  the  insulated 
wall  problem  of  Fox  and  Libby  further,  comparing  the  eigenfunction  series 
solution  with  several  other  approximate  solutions,  and  with  a  new  exact 
finite  difference  solution.  This  comparison  emphasizes  the  slow  conver¬ 
gence  of  the  eigenfunction  solution  near  the  initial  station,  and  the  remark¬ 
able  accuracy  of  a  solution  obtained  using  Lighthill's  approximate  method.  4 
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II.  FORMULATION 


In  similarity  variables  £,  r|  the  momentum  and  energy  equations  for 
unit  Prandtl  number  and  constant  density-viscosity  product  are 


f  +  ff  -  2  4  (f  f  »  -  f  t  f  )  =  0 
*1*1*1  *1*1  *1  *lS  i  *1*1 


(2.  la) 


8ii +  fgn  ‘  2eVe '  re  *  0 


(2.  lb) 


We  are  interested  in  initial  value  problems,  where  the  non-dimensional 
velocity  and  enthalpy  profiles  f^t  g  are  given  at  some  initial  streamwise 
station  £.  by  F^tq)  and  G^fq).  The  wall  boundary  conditions  are  the  usual 
ones  of  zero  velocity  and  either  constant  temperature  (problem  (2))  or  zero 
heat  transfer  (problem  (3)).  Thus 


q  =  0:  f  =  f  =0; 

*1 


q  — oe:  f^  —  1, 

i  =  if  f  =  F*(q), 


g  =  gw  for  problem  2 

g  =0  for  problem  3 
*1 


g  -  i 

g  =  G*(q) 


(2.  2) 


For  the  momentum  problem  (1),  we  perturb  about  the  Blasius  solution 


f0<*l): 


Ml  II  I  I 

f  +  f  f  =0,  f  (0)  =  f  (0)  =  0,  f  (»)  =  1 
o  o  o  ’  o'  '  o''  '  ' 


(2.  3) 


With  the  assumed  solution  as 


(2.  4) 


MM 


the  perturbation  equations  of  (2.  la)  become  for  each  f^,  k  1,2  ... 

[W  foS^fo"'-25<(o'‘,rfo"v]k  =  t'k  (2-5> 

where  F,  is  dependent  on  the  previous  perturbation  terms  f  o  <  j  <  k. 

*  J 

In  particular,  for  fj  ,  Fj  =0  and  (2.5)  is  homogenous.  In  addition,  for 
higher  perturbations,  the  homogenous  solution  can  be  used  to  find  the  par¬ 
ticular  integral.  So  we  will  focus  on  the  solution  of  the  homogenous  equation, 
which  will  lead  us  to  an  eigenvalue  problem  when  we  separate  variables. 

We  therefore  seek  solutions  of  the  homogenous  equation  asso¬ 
ciated  with  (2.  5)  in  the  form  (dropping  the  subscript  k). 


-  >_U)/2 


E  A"  (5/5i>  "  V’’1 


(2.  6) 


and  find  N  and  A  ^  must  satisfy 
n  n 


Ml  II  If  •  • 

N  +  f  N  +  A  f  N  +  (l  -  A)  f  N  =  0 
o  o  o 


(2.7) 


Since  fQ  already  satisfies  all  the  boundary  conditions  according  to  (2.  2)  and 
(2.  3)  the  boundary  conditions  on  N  are  homogeneous: 


N(0)  =  n'(0)  =  0,  n'(»)  =  0 


1 


(2.  8) 


As  discussed  in  Part  I  of  Fox  and  Libby,  (2.7)  and  (2.8)  define 
an  eigenvalue  problem,  the  eigenvalues  to  be  chosen  so  that  N'  approaches 
zero  exponentially  as  q  -*  oo.  (If  only  algebraic  decay  were  required,  there 
would  be  solutions  for  any  A.  )  To  make  the  functions  N  well-defined,  Fox 
and  Libby  imposed  the  scaling  condition 


1 


¥ 


N  (0)  =  1 


(2.  9) 


For  problems  (2)  and  (3),  f  is  taken  to  be  unperturbed  Blasius 
function  fQ(q)  in  (2.  lb)  and  the  solution  for  g  is  assumed  as 


g  =  +  g  (t.  n) 


(2.  10) 


leading  to  the  homogeneous  equation 


'  # 


g  +  f  g  2£f  gt  =  0 
Bqq  o  Bq  o  B§ 


(2.11) 


The  function  g0(Tl)  must  satisfy 


1 1  i 

g  +  f  g  =0 
Bo  o  Bo 


(2.  12a) 


i 


and  we  take  it  to  satisfy  also  the  boundary  conditions  (2.  2)  on  r|  =  0,  «: 


1=0:  S0  ■  8„ 


g  =0 
Bo 

g  1 
Bo 


for  problem  (1) 
for  problem  (2) 


(2.  12b) 


It  is  thus  completely  determined,  and  will  make  a  contribution  to  the  initial 

:|t  * 

profile  to  be  satisfied  by  g  .  However,  as  a  function  of  q,  g  satisfies  homo¬ 
genous  boundary  conditions,  since  gQ  satisfies  all  the  inhomogenous  ones, 
according  to  (2.  2)  and  (2.  12b). 

The  homogenous  equation  (2.  11)  has  separable  solutions  of  the  form 


*  ^  -\n(2'3k  (2.3) 

g  =  2J  Bn(l/U  Hn  ^ 


(2.  13) 


_  J 
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where  X^2,3^  and  H^2, 3^  satisfy 
n  n  7 

H"  +  f  h'  +  X  f  '  H  =  0  (2.  14) 

o  o 

The  boundary  conditions  on  H  for  the  constant  enthalpy  wall,  problem  (2)  , 
are  (including  a  convenient  scaling) 

Problem  (2):  H(2)(0)  =  0,  H(2)(«)  =  0;  (H(2)'(0)  =1]  (2.  15a) 

while  for  the  adiabatic  wall,  problem  (3),  they  are 

V 

Problem  (3):  H(3)'(0)  =  0,  H(3)(ao)  =  0;  [H(3)(0)  =  1  ]  (2.15b) 

These  are  two  eigenvalue  problems  with  the  same  differential  equation  but 
different  boundary  conditions  at  the  wall.  Again  the  eigenvalues  are  deter¬ 
mined  by  the  condition  that  H  -•  0  exponentially  as  q  — •  ».  These  problems 
were  discussed  by  Fox  and  Libby  in  Part  2.  ^ 

5 

Problem  (1)  has  also  been  studied  by  Stewartson,  who  showed  that 

(1)  1  /  " 

-  2  is  an  eigenvalue  with  eigenfunction  Nq  =  (q  f  -  fQ)/f q  (0).  For 

problem  (3),  Fox  and  Libby  point  out  that  =  1,  H^3^  =  f  /f '  (0), 

satisfies  (2.  14)  and  (2.  15b),  yielding  one  explicit  solution.  No  correspond¬ 
ing  result  for  either  X  or  H  is  known  for  problem  (2). 

The  standard  Sturm-Liouville  form  of  (2.  7)  is  obtained  by  recog- 

i 

nizing  that  f  is  a  solution  (though  it  does  not  satisfy  the  boundary  conditions) 

.  i  « 

so  the  introduction  of  a  new  variable  u  =  (N/fQ  )  yields 

r  *  •  Vo’2 1 » -  °  (2|6> 
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which  shows  the  weight  function  explicitly.  The  orthogonality  condition 
is  then 

/(£  ,4/£  ")  u  u  dn  =  C  (1)  6  (2.  17) 

o  o  n  m  n  mn 

with  C  ^  the  norms  for  problem  (1). 

For  problems  (2)  and  (3),  {2.  14)  can  be  put  into  standard  form  by 
1 1 

dividing  by  f  and  using  (2.  3),  to  yield 

(h7£0'Y  +  <Xf0'A0")  H  =  0  (2.18) 

and  the  orthogonality  condition  and  norms  are 

/<£  '/f  ")  H  H  dt]  =  C  (2,3)  6mn  (2.  19) 

o  o  n  m  n  mn 

The  norms  corresponding  to  the  exact  eigenfunctions  quoted  above  are  easily 
found  by  integration  to  be 

Cd1)  =  Cq3)  =  l/(2st2)  =  2.  267,  a  =  fo"(0)  =  0.4696  (2.20) 
which  both  agree  with  Fox  and  Libby's  calculated  values. 
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III.  EXTENSION  OF  KOTORYNSKI'S  APPROXIMATE 

ANALYSIS 


3 

Kotorynski  studied  these  Sturm -Liouville  problems  by  reducing 
them  to  the  form  in  which  the  coefficient  of  the  eigenvalue  is  unity.  For 
problem  (1),  this  transformation  is 


1 7  ,.2  l/4 

w  =  (f  /f  ) 

'  o  '  o  ' 


U, 


/ 


,  1/2 


t  =  /  <*  )  dr] 


(3.  1) 


which  reduces  (2.  16)  to 


d2w/dt2  +  [p-  q(t) ]  w  =  0, 


P= A 


f  f  " 

3  Vo 

q  =  4  — T~  + 

1  fit 


7 


,,2 

! _  +  I 

,3  +  4 


f' 

o 


(3.  2a) 
(3.  2b) 


For  problems  2  and  3,  the  corresponding  transformation  is 


1/4 


=  <f0'/f0"V'’  H,  t  =  J  (f0’)l/2  dr| , 


(3.  3) 


reducing  (2.  18)  to 


d2v/dt2  +  £p  -  p(t)j 


v  =  0 


i  ,  f  f  " 

1  loo 

p  =  ’  ■  p=  ‘4  7 — 


^  fo 


f 


s  f'2 

5  o  .  . 

TS  ,3  +  4  . 

o  o 


(3.  4a) 


(3.  4b) 


The  similarity  of  the  equations  is  strongly  evident,  the  functions  p,  q 
differing  only  in  the  coefficients  of  the  first  two  terms. 
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/* 


The  behavior  of  the  transformations  and  the  functions  p  and  q  is 
determined  by  the  well-known  properties  of  the  Blasius  function  fQ,  which 
are 

n  — »  0:  f  —  a  q2/2  ,  f  '  —  an  ,  f  "  -  a  ,  a  =  0.  4696  (3.  5a) 

o  o  o 


f  -  n  -  P.  f  '  -  1.  f  "  -  Ae‘(71  ‘P)  /2, 


m/  »-»; 


)  dq  =1.217 


(3.  5b) 


This  requires  that 

q  -  0:  t  -  2  (at,  3)]/2/3  (3.  5c) 

l  r  ,  1/2-,  }  f  ,  1/2. 

9  t  -  q  -  y  [l  -  (fQ  )  J  dq  -q  -J  j^l  -  (fQ  )  Jdq  (3.  5d) 

o  o 


Thus  we  find  from  (3.  2)  and  (3.  4)  that 

t  —  0:  q  — •>  7/(36t2),  p  —  -5/(36t2) 

t  q  _  p  -  (t  -  6)2/4,  6=^  ldq 


(3.6a) 

(3.6b) 


Kotorynski  proposed  a  uniform  approximation  to  p  and  q  be 
obtained  by  using  the  terms  dominant  at  t  -*0,  <*>  ,  since  the  composite 
function  would  have  the  correct  behavior  at  both  ends.  Thus  (3.  2a)  and 
(3.  4a)  are  approximated  by 


(t  -  6)2' 
4 


w  =  0, 


+ 


(3.  7a,  b) 
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Actually,  Kotorynski  did  not  include  the  6  term  in  these  equations. 
Without  it,  they  are  solvable,  as  he  showed.  With  it,  they  are  apparently 
not  solvable  in  terms  of  known  functions.  However,  the  application  of  the 
WKB  method  for  finding  the  eigenvalues  is  possible  even  when  6  is  retained, 
and  the  6  terms  greatly  improve  the  accuracy  of  the  resulting  eigenvalues, 
compared  to  those  found  by  Kotorynski,  as  we  shall  see  in  Section  V. 

A  comparison  between  the  exact  functions  q  and  p  of  (3.  2b)  and  (3.  4b) 
and  the  approximations  given  in  (3.  7a,  b)  is  shown  in  Fig.  1.  They  are  plotted 
against  q  ,  which  is  related  to  t  by  (3.  1).  The  approximations  are  shown  for 

5 

both  6=0  and  6  =  0.  425,  which  is  the  value  given  in  Stewartson.  For  q, 
there  is  a  small  difference  between  the  exact  function  and  the  approximation 
with  6  \  0  for  0.  6  <  q  <  3.  2,  but  they  agree  quite  well  elsewhere.  If  6  =  0, 
the  approximation  is  substantially  in  error  for  q  >  3.  For  p,  there  is  very 
little  error  in  the  approximation  when  6  ^  0  to  the  scale  of  Fig.  1;  if  6  =  0 
there  are  again  substantial  errors  for  r|  >  3.  It  is  clear  that  using  the 
correct  asymptotic  representation  of  fQ  makes  a  large  improvement  in  the 
accuracy  of  approximation  of  p  and  q. 

Before  the  approximate  equations  can  be  solved,  the  boundary 
conditions  must  be  stated  in  terms  of  w  and  v.  If  one  uses  the  boundary 
conditions  (2.  8)  and  (2.  9),  traces  through  the  transformation  to  u,  (2.  16), 
and  to  w,  t  (3.  1),  and  uses  (3.  5c),  one  finds 

t  -  0:  w  —  (3t/2)7/6/2  a1/3  (3.  8a) 
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Fox  and  Libby*  have  shown  that  the  asymptotic  form  of  N'  is 
N'  -  (tl-P)X_1  exp  £-(r|  -  (3)2/2  J 
When  this  is  transformed  to  w,  t  using  (3.  5b,  d),  we  find 


t  -*  oo  ;  w  ~  (t-6)X  "*  exp 


[-<t-6)2/4] 


(2)  (3) 

Similarly,  the  boundary  conditions  (2.  15a)  and  (2.  15b)  for  v  and  v 
yield,  respectively,  for  problems  2  and  3: 


t  —  0:  t-1/6  v(2)  —  0,  (3  at/2)1/6  ^ 


(2)  _  ..(2) 


dt 


a  v’  ~/4  1 

(3  a  t/2)5/6 


r'3»  -  (3  t/2  a2)1/6,  (3  at/2)1/6 


dt 


(3  t/2) 


273 


t  —  :  v  ~  (t-6)X  1  exp 


[-(t-6)2/4] 


We  may  now  proceed  to  the  solution  of  (3.  7)  for  6=0,  A  change 

l/2 

of  variable  t  =  2  '  x  permits  them  both  to  be  written 


d2  y/dx2  +  (  2(ji  -  x2  +  (4  1  -  q2)/x2]  =  0 
Problem  1:  a  =  _+  2/3;  Problem  2,  3:  a  =  +  l/3 
and  a  solution  of  this  equation  is  given  in  Ref.  6,  p.  781,  Eq.  (22.  6.  18) 

as  y 

-x  /2  q  +l/2  T  (q  )  ,  2. 
y  ~  e  '  x  '  Ln'  (x  ), 

2p  =  2A-l  =  4n  +  2a  +  2,  n=0,  1,  2,  ... 


where  L  (q  )  (x2 
n 

Ref.  6,  p.  775, 


)  is  the  generalized  Laguerre  polynomial  defined  in 

2 

Eq.  (22.  3.  9),  as  an  n-th  order  polynomial  in  x  .  In 


(3.  8b) 


* 


(3.  9a) 


(3.  9b) 


(3.  9c) 


(3.  10) 


(3.  11) 
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terms  of  t,  we  may  take  as  our  solutions* 


2 

L  ,  v  «>.  v  <3>1  (d  <*>.  D  <2>,  D  (3))  e''  74  t“  +  1/2L  <“»(,2/2) 

with  the  to  be  determined  by  the  boundary  conditions. 

First  we  notice  that  this  solution  satish  s  the  boundary  condition 
as  t  -»«•  ;  in  fact  (3.  12)  becomes  proportional  to  e  1  ^  t2n  *  Q  *  *^2 
which,  from  (3.  11),  agrees  exactly  with  (3.  8b)  and  (3.  9c).  Second,  we 
find  the  behavior  of  w,  v  and  their  derivatives  as  t  -*  0  from  Ref.  6, 
p.  775,  Eq.  (22.  3.  9): 


(3.  12) 


t  —  0:  w,  v  ~  ta  +  1//2  [T  (n  +  a  +  l)/n!  T  (  a  +  1)  +  O  (t2)]  (3.  13a) 

w’,  v’  ~  (a  +  1/2)  ta_1//2  [T  (n  +  a  +  l)/n!  T  (  a  +  1)  +  O  (t2)]  (3.  13b) 


Third,  the  eigenvalues  are  given  in  terms  of  a  by  (3.  11).  For  each 

problem,  there  are  two  possible  values  of  a,  Eq.  (3.  10),  only  one  of 

which  will  satisfy  the  boundary  conditions  at  t  =  0. 

Kotorynski  did  not  actually  exhibit  the  resulting  approximate 

eigenfunctions  which  satisfy  the  boundary  conditions  by  evaluating  the 

constants  D  ,  nor  did  he  find  the  norms  C  .  However,  the  expression 
n  n 

for  the  norms  is  easily  found  from  the  orthogonality  properties  of  the 
Ln^a).  Comparison  of  the  definitions  (2.  17)  and  (2.  19)  with  the  trans¬ 
formations  (3.  1)  and  3.  3)  show  that  we  need  the  integral 

ij  e-'2/2  .2o+1  rLn'“>  («2/2> 

0  *- 


*The  solution  given  in  (3.  12)  can  be  compared  with  Kotorynski' s  Eq.  (3.  5). 
Our  a  is  his  2a,  but  his  argument  of  L  does  not  have  the  factor  l/2.  This 
represents  either  a  different  definition  of  L  (which  he  does  not  give)  or 
an  error  in  his  formula. 
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But  Eq.  (22.  2.  12)  on  p.  775  of  Ref.  6  shows  that  with  z  =  t  /2, 


2  "  Q  I 


=  J  e  z  za  ^Ln^a  ^  (z)]  dz  =  r  (n  +  a  +  l)/n! 


Thus  the  norms  are  related  to  the  coefficients  D  by 

n 

C  =  2a  D  2  r  (n  +  a  +  1  )/n ! 
n  n 


(3.  14) 


The  approximate  solutions  can  now  be  found  by  choosing  to 
satisfy  the  boundary  conditions  (3.  8a),  (3.  9a,  b)  as  t  -»  0,  using  the 
asymptotic  forms  from  (3.  13a,b).  In  the  process  of  satisfying  the  boundary 
conditions,  one  finds  which  one  of  the  two  a  values  for  each  problem, 

Eq.  (3.  10),  is  the  correct  one;  this  determines  the  approximate  value 
of  the  eigenvalues  from  Eq.  (3.  11).  The  norms  follow  from  Dn  by 
(3.  14). 


The  results  of  applying  this  process  can  be  expressed 


compactly; 


„(1)  n<2)  n<3) 

a  ,  a  ,  a 


2/3,  1/3,  -1/3 


l(i)  +  1/2  =  A(l)  =  2n  +  a(l)  +  3/2 


(3.  15a) 


(3.  15b) 


\  m  (2)  (3)1  /  3  ^  1  +l/z  nir(q(i)+  1) 

jD<"  K  •  Dn'  j--(2a2)  r <«'*'+  n+1)  2  '  *' 


(3.  16) 


(i)  _  n!  T(au4  1) 


r(a(l)  +  n  +  1) 


C  ^  -  n  c(i) 

o  „  ,  „  (i)  n-1 

n  +  a 


jc(1)  ,  C  (2), 

1  o  o’ 


Co(3>  =  2 


(2a2) 


2a(i)+l 


T(a(l)+  1) 


1.  52,  6.  05,  2.  04 


(3.  17a) 


{*• a2'  '} 


(3.  17b) 
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This  completes  the  solution  of  the  three  problems  when  q,  p  are  approximated 
by  the  sum  of  their  asymptotic  expressions  as  t  -»  0  and  00  ,  and  the  constant 
6  in  the  latter  asymptote  is  taken  to  be  zero. 
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IV.  CONVERGENCE  OF  HEAT  TRANSFER  SOLUTIONS 


We  will  now  consider  briefly  the  solutions  problems  2  and  3, 
concentrating  on  the  convergence  of  the  series  solutions.  For  this  purpose, 
the  approximate  eigenvalues  and  normals  obtained  in  the  previous  section 
will  be  used.  Although  they  are  not  particularly  accurate  for  moderate  n, 
they  become  accurate  for  large  n,  and  so  are  useful  for  convergence 
studies. 

In  problem  2,  Fox  and  Libby ^  consider  a  flat  plate  which  is  held 

at  a  fixed  enthalpy  gw>  up  to  and  then  changed  to  gWQ.  A  function  gQ 

which  satisfies  (2.  12)  is  gQ  =  gwQ  +  0"gWQ)  f  •  Then  g*  is  given  by  (2.  13) 

(2) 

and  (2.  14)  with  boundary  conditions  (2.  15a).  The  coefficients  are 

determined  by  the  initial  condition  on  g,  which  is  that  up  to  £^,g  is  the 

I 

Crocco  integra  1  8  =  Swl  +  (l-gwi)  f  .  The  initial  condition  on  g*  then 

(2) 

determines  as  an  expansion  coefficient: 

(g  -  -  g  )  (1-f  ')  =  £  b(2)  h(2) 

'Bwl  "wo  '  o  4-*  n  n 

o 

Using  the  orthogonality  condition  (2.  19)  as  well  as  the  Sturm- Liouville 
form  of  the  equation  for  H,  (2.  18),  the  integral  can  be  evaluated  to  give 

B(2)  =  (g  -g  )/fx(2)  a  C(2))  a  =  f  "(0)  =  0.  4696 

n  &wl  6wo  '  V  n  n  /  o 


The  heat  transfer  at  the  wall  is  found  using  this  in  (2.  13)  and  differentiating 


(2.  10): 


g  =  ( 1  - g  )a  +  (g  ,-g  )a 
Br)W  Bwo  ewl  6wo 


1  i(^] 


(4.  1) 
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This  agre  with  Fox  and  Libby,  Part  2,  Eqs.  (2.  26b)  and  (2.27),  who 
take  an  initially  adiabatic  wall  so  g^  =  1. 

In  problem  3.  the  wall  is  still  held  at  g^  up  to  £.,  but  then  changes 
to  an  adiabatic  wall,  g  =0.  In  this  case  g  =1.  and  the  initial  condition 
on  g  is  the  same  Crocco  integral  as  above.  Therefore,  the  initial  condition 
for  this  problem,  using  again  the  solution  (2.  13)  for  g-; .  is 


-(1-*wl> 


(1 


f  ) 

o 


E 

o 


(3)  (3) 

B  H 
n  n 


The  coefficient  is  found  the  same  way  as  before,  although  the  result  differs, 
since  the  H  satisfies  different  boundary  conditions. 


B(3)  =  -Cl -B  ,)  /  (A  (3)  C(3)  ) 

n  *wr  '  n  n 


Using  this  in  (2.  13),  the  wall  temperature  is  found  from  (2.  10)  as 

oo 


*  E  <Xn’  Ci31' 


(^/4i)*Xn  /2 


(4.  2) 


This  agrees  with  Fox  and  Libby,  Part  2,  Eqs.  (2.  30b)  and  (2.  31),  when 

a  misprinted  subscript  o  is  changed  to  1  in  the  latter  equation. 

We  now  consider  the  convergence  of  the  sums  in  (4.  1)  and  (4.  2). 

using  (3.  15)  and  (3.  17)  for  A  and  C  .  The  ratio  of  the  n  +  1  to  the  n  term, 

n  n 

for  large  n,  becomes 

r  =  [l  +  (a(l)-l)/n]  £./£ 

For  £>£.,  this  converges  by  the  ratio  test,  although  the  convergence  is 
very  slow  near  £..  For  £  =  £^,  at  the  initial  plane,  the  series  only 
converges  if  -  1  <  -1,  i.  e.  ,  a  ^  <  0.  Thus  the  series  (4.  2)  for 
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/  Q  \ 

problem  3  converges  at  §  =  £.,  since  a  =  l/3,  but  the  series  (4.  1)  for 
problem  2  diverges,  since  a  -  l/3.  This  is  in  accord  with  the  indications 
of  the  numerical  results  in  Fox  and  Libby,  Part  2,  where  Fig.  4  shows 

rising  steeply  near  the  initial  plane,  while  Fig.  5  has  a  finite  value  for 
g  there. 

Notice  that  if  one  differentiates  (4.  2)  with  respect  to  £  ,  one  brings 

down  a  X  ,  cancelling  that  term  in  the  denominator.  Then  the  ratio  of 
n  ° 


terms  for  large  n  is 


r  =  (1  +  a(l)/n)  %/%. 


which  still  yields  convergence  for  £  >  £^,  but  now  shows  that  problem  3 

/  O  \ 

has  a  divergent  streamwise  derivative  of  g  ,  since  a  =  -l/3  >  -1. 

This  would  indicate  that  the  curves  of  Fig.  5  in  Fox  and  Libby,  Part  2, 
should  have  an  infinite  slope  at  the  initial  plane  a  property  that  their  solution 
does  not  exhibit  in  the  Figure.  This  is  likely  caused  by  the  use  of  only  10 
terms  of  a  slowly  -  converging  series.  Section  VI  has  a  fuller  discussion 
of  problem  3,  where  it  is  shown  that  the  next  10  terms  make  a  substantial 


contribution,  and  that  the  slope  at  the  initial  plane  is  indeed  very  steep 
and  probably  infinite. 
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V.  ESTIMATION  OF  EIGENVALUES  BY  THE  WKB  METHOD 


The  eigenvalues  of  (3.  7a,  b)  can  be  estimated  by  a  technique  well 
known  in  quantum  mechanics  for  handling  the  Schrttdinger  equation.  The 
advantage  of  this  over  Kotorynski's  approach  is  that  the  correct  behavior 
for  large  t,  which  retains  the  constant  6,  can  be  used.  This  leads  to  a 
much  more  accurate  estimate,  but  does  not  easily  give  approximate 
eigenfunctions.  However,  once  accurate  eigenvalues  are  known,  the 
numerical  integration  of  the  exact  equations  (2.  7)  or  (2.  14)  presents 
little  difficulty.  The  difficulties  encountered  by  Fox  and  Libby  are 
primarily  associated  with  the  necessity  of  finding  the  eigenvalues  by 
trial  and  error. 

The  idea  of  the  WKB  method  is  to  solve  (3.  7a,  b)  in  three  regions. 

-2  2 

One  is  where  t  is  very  large,  so  the  t  term  is  negligible  and  (t-6)  /4 >p; 
in  this  region  there  are  exponential-type  solutions.  A  second  region  is 
where  t  is  still  large,  so  t  is  neglected,  but  (t-6)  /4  <  p  so  the  solutions 
have  sinusoidal  character.  These  two  solutions  are  matched  at  the  turning 
point  t  where 

(tQ-6)2/f=H  (5.D 

so  as  to  eliminate  the  increasing  exponential  term,  and  make  the  solution 

continuous.  (This  matching  procedure  is  described  in  detail  in  Ref.  6, 

for  example.  )  To  satisfy  the  boundary  conditions  at  t  =  0,  a  third  solution 

“2  2 

is  found  in  the  region  where  the  t  term  dominates  (t-6)  /4.  Once  the 
boundary  conditions  are  satisfied  with  this  solution,  its  asymptotic 
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expansion  is  joined  to  the  solution  in  the  second  region;  this  joining  gives 

a  condition  for  determining  the  eigenvalues. 

This  procedure  can  be  carried  out  in  general  for  (3.  7a,  b)  and  the 

appropriate  boundary  conditions  (3.  8a)  or  (3.  9a,  b).  However  the  result 

is  already  available  in  a  quadrature  formula  in  Ref.  7,  as  pointed  out  by 

Kotorynski,  who  used  it  with  6  =  0  to  verify  his  approximate  eigenvalues 

7 

(3.15).  Titchmarsh  gives,  on  page  151,  Eq.  (7.76),  a  formula  which 
for  (3.  7a,  b)  becomes 

1/2 

dt  =  +  n  +  0(— ),  n=  0,1,2 _ 

c.  n 


where  tQ  is  defined  in  (5.  1),  By  integrating  the  left  side,  we  obtain  an 
equation  for  p^. 

Carrying  at  the  integration  by  a  trignometric  substitution  we  find 


This  shows  pn  =  O(n),  so  we  may  keep  terms  in 


The  result  is 


x  26  1/2 

H-  +  —  H  ' 


1 2  7r  ^n 


1/2  _ 


+  l  +  2n  +  O  (— ) 
n 


1  /2 

This  is  a  cubic  in  p^  '  ,  but  the  last  term  on  the  left  is  quite  small,  because 

5 

Stewartson  has  given 
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Therefore,  we  will  ignore  this  term,  and  solve  the  quadratic  to  get 
jin^2  =  Un-  1/2)1/2  =  -  6/j r  +  [(  6/tt)2  +  a  +  1  +  2n]*  ^ ,  n  =  0,  1,2,...  (5.  3) 

For  6=0,  (5.  3)  agrees  exactly  with  (3.  15),  which  was  the  result 
obtained  by  Kotorynski.  The  improvement  obtained  by  including  the 
6  term  is  remarkable,  as  shown  in  Table  I,  where  the  exact  results  of 
Fox  and  Libby  are  compared  with  those  of  Kotorynski,  Eq.  (3.  15),  the 
present  results,  Eq.  (5.  3),  and  some  additional  exact  results  for  problem 
3,  obtained  by  the  present  author.  Even  at  the  lowest  eigenvalues,  the 
error  is  only  7.  4%  for  problem  1,  1.  1%  for  problem  2,  and  2.  1%  for 
problem  3.  This  accuracy  is  a  pleasant  surprise  since  the  WKB  method 
is  designed  for  the  higher  eigenvalues.  The  errors  are  less  than  1% 
beginning  with  the  fourth  eigenvalue  for  problem  1,  and  the  second  for 
problems  2  and  3;  they  are  less  than  0.  1%  after  the  thirteenth  eigenvalue 
for  problem  1  and  the  fourth  for  problems  2  and  3.  It  is  clear  that  the  use 
of  (5.  3)  will  enable  the  numerical  solution  of  (2.  7)  or  (2.  14)  to  be  performed 
in  one  forward  integration  from  rj  =  0. 

If  Eqs.  (3.  7a,  b)  could  be  solved  exactly  with  the  6  terms  included, 
presumably  very  accurate  norms  could  also  be  found,  but  the  present 
author  has  been  unable  to  find  such  a  solution  in  terms  of  known  functions. 

The  norms  obtained  by  the  Kotorynski  approximation,  (3.  17),  are  shown  for 
problem  3  in  Table  II,  together  with  exact  values  calculated  by  the  present 
author.  The  errors  vary  from  10%  at  n  =  0  to  1%  at  n  =  19.  (The  first 
10  exact  values  of  C  in  Table  II  differ  in  the  second  decimal  place  from 
those  given  on  p.  438  of  Fox  and  Libby,  ^  Part  2.  They  are  believed  to  be 
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TABLE  I.  1 

Eigenvalues  for  Problem  1,  n  =  2 


Xn(exact)2 

Xn(Eq.  3.  15)3 

Xn(Eq.  5.3) 

n 

Xn(exact)2 

Xn(Eq.3.  15)3 

Xn(Eq.5.  3) 

2 

2.  167 

1. 852 

10 

20.  984 

22.  167 

20. 943 

3.  774 

4.  167 

3.  684 

n 

22.924 

24.  1  67 

22. 886 

5.  627 

6.  167 

5.  558 

12 

24. 862 

26.  167 

24. 832 

7.  517 

8.  167 

7.  453 

13 

26.  865 

28.  167 

26.  780 

9.408 

10.  167 

9.  361 

14 

28.  755 

30.  167 

28.729 

11. 328 

12.  167 

11.  278 

15 

30. 711 

32.  167 

30.  680 

13. 252 

14.  167 

13.  202 

1  6 

32.  667 

34.  167 

32. 633 

15. 171 

1  6.  1  67 

15. 132 

17 

34.  620 

36.  1  67 

34. 587 

17.099 

18.  167 

17. 066 

18 

36.  567 

38.  167 

36.  542 

19.040 

20.  167 

19. 003 

19 

38.  520 

40. 167 

38. 499 

TABLE  I.  2 


Eigenvalues  for  Problem  2,  a  =  1/3 


n 

X^exact)1 

An(Eq.  3.  15)3 

Xn(Kq.  5.3) 

n 

Xn(  exact)1 

Xn(Eq.  3.  15)3 

Xn(Eq.  5.3) 

0 

1.  572 

1.  833 

1.  555 

5 

10.  96 

11. 833 

10. 958 

1 

3.  385 

3.  833 

3.  375 

6 

12.  88 

13. 833 

12.  881 

m 

5.25 

5.  833 

5.  244 

7 

14.  81 

15.  833 

14. 810 

| 

7.  14 

7.  833 

7.  136 

8 

16.  74 

17. 833 

16.  743 

B 

9.  05 

9.  833 

9.  042 

9 

18.  68 

19.  833 

18. 680 

TABLE  I.  3 
Eigenvalues  for  Problem  3,  0  =  1/3 


n 

X  (exact)3 
n 

An(Eq.  3.  15)3 

An(tq.  5.3) 

n 

XR(  Exact)3 

XJEq.  3.  15)3 

An(Eq.  5.3) 

0 

i 

1.  167 

0.979 

10 

19. 978 

21.  167 

19. 973 

1 

2.774 

3.  167 

2.  760 

11 

21.920 

23.  167 

21. 915 

4.  629 

5.  167 

4.  618 

12 

23.  864 

25.  167 

23.  859 

■ 

6.  513 

7.  167 

6.  504 

13 

25.  811 

27.  1  67 

25. 806 

8.  414 

9.  167 

8.406 

14 

27.  759 

29.  1  67 

27. 754 

5 

10. 326 

11.167 

10. 319 

15 

29.  709 

31.  167 

29.  704 

6 

1  2.  247 

13.  167 

1  2.  240 

16 

31.  661 

33.  167 

31.  656 

7 

14.  173 

15.  167 

14. 1 66 

17 

33.  614 

35.  167 

33. 610 

8 

16.  104 

17.  167 

16.  098 

18 

35.  569 

37.  167 

35.  565 

.  18.040 

19.  167 

18. 034 

19 

37.  525 

39.  1  67 

37. 520 

£ 

These  exact  results  have  been  calculated  by  the  present  author.  The  first  10  agree 
with  Fox  and  Libby,  Part  2,  to  within  0.  01 


TABLE  II 

Norms  for  Problem  3,  a  =  -1/3 


n 

Cn(exact) 

Cn(Eqs.  3.  1  7)3 

n 

Cn(  exact) 

Cn(Eqs.  3.  1  7)3 

0 

2.  268 

2.  035 

10 

6.  077 

6.  00 

1 

3.  184 

3.045 

1 1 

6.  262 

6.  19 

2 

3.  772 

3.  650 

12 

6.  438 

6.  36 

3 

4.  219 

4.  11 

13 

6.  604 

6.  53 

4 

4.  586 

4.49 

14 

6.  763 

6.  69 

5 

4.903 

4.  81 

15 

6.915 

6.  84 

6 

5.  182 

5.09 

16 

7.059 

6.  99 

7 

5.  435 

5.  35 

17 

7.  198 

7.  12 

8 

5.  665 

5.  58 

18 

7.  333 

7.  26 

9 

5.  878 

5.  80 

19 

7.  461 

7.  39 
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* 


g 

more  accurate,  since  the  method  of  Nachtsheim  and  Swigert  was  used 
to  obtain  them,  and  the  criterion  for  stopping  the  integration  was  the 
constancy  of  itself  as  r)  —  ) 
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VI.  THE  INSULATED  WALL  PROBLEM 


Problem  3,  a  flat  plate  with  constant  wall  condition  g^  up  to 

£.  and  insulated  (e  =  0)  for  £  >  £.,  has  had  some  previous  attention. 

*1  sr|w  *  n  r 

9  10 

In  particular,  Durgin  ’  made  a  thorough  experimental  and  theoretical 

study,  reviewing  the  previous  solutions,  all  of  which  were  approximate. 

He  concluded  that  the  laminar  experimental  results  for  wall  temperature 

distribution  in  the  insulated  region  were  in  very  good  agreement  with  the 

4 

theory  based  on  Lighthill's  approximate  method,  thus  giving  it  credence 
as  the  best  approximation  to  the  exact  solution.  The  solution  of  Fox  and 
Libby,  ^  given  in  (4.  2),  is  an  exact  solution  expressed  in  series  form  and 
may  be  compared  with  the  Lighthill- Durgin  solution. 

Strangely  enough,  no  numerical  solution  to  the  partial  differential 
equation  of  the  problem  seems  to  have  been  published  up  to  this  time, 
although  such  a  solution  would  represent  another  exact  solution  to  serve 
as  a  standard  of  comparison  not  only  for  the  Lighthill  -  Durgin  approximation, 
but  for  the  rate  of  convergence  of  the  Fox- Libby  series. 

Such  a  numerical  solution  has  been  obtained  as  a  by-product  of  a 
computer  program  developed  for  another  problem.  ^  ^  It  solves  the  energy 
equation  derived  from  (2.  2)  when  f  =  fQ(r|): 

g  +  f  g  -  2£  V  g,  =  0  (6.  1 ) 

hr|  r)  obr|  *  o  *£  '  ' 
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The  boundary  and  initial  conditions  are: 


S  >  g  (n  =  0)  =  0,  g(n  — 00  )  —  1 


^  ^i*  ®i  "  ®wl  *  <1-8wl)  *o 


(6.  2a) 
(6.  2b) 


The  parameter  g^  is  eliminated  by  the  transformation 

G=(g-D/(gwl-  1) 

and  the  streamwise  coordinate  is  simplified  by  introducing  a  new 


(6.  3) 


coordinate 


s  =  In  (S/Sj) 


1/2 


z  >  0  for  £  >  £. 


(6.  4) 


Then  (6.  1)  and  (6.  2)  become 


G  +  f  G  -  f  G  =  0 
tin  o  T|  o  z 


z  >  0:  G^  (t)  =  0)  =  0,  G  (r|  -•*>)  =  0 


z  =  0:  G.  =  1  -  f 
1  o 


(6.  5) 

(6.  6a) 

(6.  6b) 


This  is  a  problem  with  no  parameters.  The  new  streamwise  coordinate 
z  seems  to  be  the  natural  one  for  initial  value  problems  (where  £.  *  0) 
because  it  takes  the  operator  2  £  9/9  £  into  9  /d  z.  It  is  rather  surpris¬ 
ing  that  it  does  not  seem  to  have  had  any  widespread  use  in  numerical 
solutions  of  boundary  layer  problems. 

The  solution  to  (6.  5),  (6.  6)  was  obtained  by  use  of  a  Crank- 

12 

Nicolson  implicit  finite  difference  scheme.  The  r|  derivatives  are 
averages  of  central  difference  formulas  at  the  known  station  z  and  the 
new  station  z  +  Az,  while  the  z  derivative  is  a  forward  difference.  This 
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leads  at  each  z  to  a  set  of  linear  inhomogenous  algebraic  equations,  each  of 

which  contains  the  unknown  at  only  3  successive  q  points.  The  matrix  of 

coefficients  is  in  tri-diagonal  form,  and  the  equations  are  easily  solved 

by  successive  elimination,  starting  at  q  =  0,  moving  to  the  edge  of  the 

boundary  layer,  and  then  sweeping  back  toward  the  wall  again.  The 

12 

algorithm  is  easily  formalized. 

The  boundary  condition  G  =  0  at  q  -»  00  is  satisfied  at  a  finite 
value  q  ,  chosen  so  that  f  is  very  near  its  asymptotic  value  of  unity.  We 
used  qe  =  4.  8  as  the  edge.  The  wall  condition  is  on  which  presents 
some  complications,  since  the  central  difference  at  the  wall  involves  a 
point  inside,  while  that  one  step  off  the  wall  involves  the  unknown  G  . 

But  the  differential  equation  (6.  5)  at  the  wall,  and  the  boundary  condition 
(6.  6a)  together  tell  us  that 


G 


qqw 


=  0 


Therefore,  a  Taylor  series  expansion  from  the  wall  q  =  0  to  the  first 
station  q  =  A  q  shows  that 

Gw  =  G(Aq)  +  0(Aq2) 

To  the  order  of  accuracy  of  the  central  difference  formulas  then.  G^  =  G(Aq), 

and  this  relation  makes  the  algebraic  solution  of  the  difference  equations 

determinate,  if  we  start  with  the  equation  centered  at  q  =  Aq  and  end  with 

the  one  at  q  =  q  -  Aq,  taking  G(q  )  =  0. 

By  this  scheme,  solutions  have  been  obtained  with  q^  =  4.  8,  Aq  = 

-3  -4  -5 

0.  05,  for  Az  =  5x10  ,  5x10  ,  and  5x10  .  The  coefficient  functions,  f  , 

o 

f'  were  generated  first,  by  solving  (2.  3)  with  f  "(0)  =  a  =  0.  4696  with  the 
same  A  q.  The  results  for  G  are  shown  in  Fig.  2  over  the  range  O^zSTO.  13, 
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which  corresponds  to  1<£/£^<1.  3.  This  short  distance  near  the  initial 

plane  has  been  singled  out  to  study,  because  the  behavior  near  the  initial 

plane  is  the  most  stringent  test  of  any  approximate  solution. 

The  numerical  solutions  seem  to  approach  =  1  smoothly,  but 

with  a  very  steep  slope,  which  can  well  be  interpreted  as  infinite  at  z  =  0. 

The  Fox-Libby  series,  (4.  2),  is  also  shown  through  10  terms  by 

the  X  symbols.  They  show  good  agreement  with  the  numerical  solutions 

(within  10%)  for  z  >  0.  09,  but  fall  below  it  for  smaller  z,  reaching  only 

=  0.  7426  at  z  =  0.  (This  solution  was  calculated  using  the  exact 

eigenvalues  given  in  Table  I.  3,  and  the  exact  norms  shown  in  Table  II.  ) 

Extending  the  series  to  20  terms  yields  the  results  shown  on  Fig.  2  by  the 

solid  diamonds.  The  solution  is  improved,  but  still  only  reaches  G  =  0.  8 

w 

at  z  =  0.  Clearly,  many  terms  are  needed  to  obtain  accuracy  near  the 
initial  plane  because  of  the  slow  convergence  of  the  series  (4.  2),  although 
we  proved  in  Section  V  that  it  did  converge  and  could  presumably  yield 
Gw  =  1  at  z  =  0  eventually. 

The  solution  to  this  problem  based  on  Lighthill's  method  was  worked 

9 

out  by  Durgin.  The  velocity  profile  is  approximated  by  a  linear  function 
f'  =  ar|,  VonMises  variables  are  used,  and  the  equation  is  solved  by  a 
Laplace  Transform  in  the  streamwise  coordinate.  One  may  then  invert 

4 

back  to  the  physical  plane  in  two  ways.  Lighthill  finds  heat  transfer  as 
an  integral  over  the  wall  temperature  distribution.  For  the  present 
problem  we  want  the  wall  temperature  as  an  integral  over  the  heat 
transfer.  Into  this  integral  we  substitute  the  shear,  and  the  heat 

transfer  of  the  similar  boundary  layer  over  the  front  of  the  plate  0<^<^. 
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For  this  latter  heat  transfer  we  use  the  value  obtained  from  the  Lighthill 
approximation  for  a  constant  temperature  wall,  to  be  consistent  with  the 
approximation  being  used  for  The  result  of  integration  is  then, 

for  $>$., 


./e)374 

du 

1/3  .2/3 

u  '  (1 -u)  ' 

9 

which  was  given  by  Durgin. 

While  this  is  an  incomplete  Beta  function,  its  tabulation  can  be 
easily  accomplished  by  using  Simpson's  rule  on  the  integral,  after  re¬ 
moving  the  singularity  of  the  integrand.  To  accomplish  this  we  let 


(6.  7) 


=  (cos  9  ) 


-8/3 


and  find 


2  sin  7r/ 3 
7 r 


(6.  8) 


The  result  of  such  a  calculation  is  shown  in  Fig.  2  as  the  +  symbols. 
They  show  remarkable  agreement  with  the  numerical  solution  right  up 
to  z  =  0,  including  the  rapid  increase  in  slope.  From  (6.  7)  the  slope 
near  z  =  0  is 


dG 


_ w 

dz 


(3/2)*^  sin  7r/3 


7T  Z 


2/3 


0.  316 

273 

z  ' 


(6.  9) 


which  shows  the  strength  of  the  singularity, 
as  i/e.  -*  00  ,  we  find  from  (6.  7)  that 


3  sin  ?r/3 
2  7T 


t  -1/2 

If-)  =  0.  414 

^•5 


-1/2 

) 
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which  may  be  compared  with  the  corresponding  result  from  the  Fox-Libby 
series,  which  is  just  the  first  term,  i.  e. 

Gw  -  (!/|.)-l/2/2.  268  =  0.  440  (S/l.)’1/2  (6.  10) 

Since  the  latter  result  is  presumably  exact,  we  see  the  Lighthill- 
Durgin  solution  is  some  6%  low  for  large  £.  So  its  accuracy  is  likely 
to  be  high  over  the  whole  range. 

Durgin  already  found  (6,  7)  to  be  in  good  agreement  with  his 
experimental  results,  and  now  we  have  shown  it  to  be  also  in  good  agree¬ 
ment  with  our  exact  numerical  solution.  This  leaves  little  doubt  that  the 
simple  integral  (6.  7)  is  the  most  useful  available  solution  to  the  present 
problem.  One  may  speculate  on  why  a  solution  developed  for  large 
Prandtl  number  should  be  so  accurate  for  unity  Prandtl  number.  We 
already  know  that  the  Lighthill  approximation  for  heat  transfer  predicts 
nearly  the  correct  Prandtl  number  dependence,  the  l/3  power,  although 
a  somewhat  erroneous  constant.  But  the  relation  leading  to  (6.  7)  involves 
the  heat  transfer  divided  by  the  l/3  power  of  the  Prandtl  number,  so  in 
fact,  (6.  7)  is  obtained  as  independent  of  the  Prandtl  number,  and  holds  for 
any  value  of  that  parameter.  It  is  likely  to  give  a  good  answer  when  com¬ 
pared  with  numerical  solutions  obtained  for  nonunity  Prandtl  numbers. 

(In  fact,  Durgin' s  experiments  were  for  air,  which  has  a  Prandtl  number 
about  0.  72.  )  One  must  only  be  careful  to  define  G  to  be  (T  -T  )/(T  ,-T  ) 
where  is  the  adiabatic  wall  temperature,  so  the  correct  limit  is 

approached  for  downstream.  In  the  present  case  of  unit  Prandtl  number, 

Taw  is  the  same  as  the  free  stream  stagnation  temperature,  so  our 

definition  of  G  ,  (6.  3),  is  consistent  with  this, 
w 


-35- 


As  an  example  of  integral-type  solutions,  we  may  quote  that  of 
1  3 

Libby  and  Morduchow,  which  Fox  and  Libby  also  plot  in  their  Part  2, 

9 

Fig.  5,  and  Durgin  quotes.  For  Prandtl  number  unity  and  constant 
density-viscosity  product,  this  result  is 

gw  =  o.476  (i/e.)'1/2  +  o.  524  (e/e.r7, 63 

This  is  also  plotted  as  the  triangle  symbols  in  Fig.  2.  It  clearly  has 
the  wrong  behavior  near  z  =  0,  and  is  8%  high  as  £  -*  °°  according  to 
(6.  10).  In  the  intermediate  range,  starting  about  =  2.  5,  it  is  in 

rather  good  agreement  with  10  terms  of  (4.  2),  better  than  Fox  and  Libby 
show  in  their  figure,  where  they  seem  to  have  an  error  in  the  Libby- 
Morduchow  curve.  But  in  any  event,  it  is  less  accurate  than  (6.  7),  and 
hardly  any  more  convenient. 
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VII.  CONCLUSIONS 


By  pursuing  further  Kotorynski' s  uniform  approximation  it  has  been 

» 

possible  to  obtain  explicit  expressions  for  the  eigenfunctions  and  their  norms 
in  three  boundary  layer  eigenvalue  problems.  The  convergence  of  the  Fox- 
Libby  eigenfunction  series  solution  to  two  heat  transfer  problems  has  then 
been  investigated,  using  these  norms  and  the  eigenvalues  obtained  by 
Kotorynski.  The  series  both  converge  away  from  the  initial  plane  with  a 
convergence  factor  but  at  the  initial  plane,  the  series  for  a  discontin¬ 

uous  wall  temperature  diverges.  The  series  for  an  insulated  wall  following 
a  constant  temperature  wall  converges  at  the  initial  plane,  but  its  stream- 
wise  derivative  diverges  there. 

We  have  shown  that  the  approximate  eigenvalues  obtained  by  Kotorynski 
)  can  be  greatly  improved  by  using  the  WK  B  method.  This  improvement  is 

obtained  by  using  a  more  accurate  approximation  to  the  Blasius  fund  on  for 
large  r|,  i.  e.  ,  using  fQ  —  r|-P  instead  of  f  —  r|.  While  the  Kotorynski  type 
of  explicit  solution  cannot  be  obtained  if  (3  is  kept,  the  WKB  method  for  finding 
the  eigenvalues  can  be  easily  applied.  The  resulting  eigenvalues  are  iound 
as  the  solution  of  a  quadratic  equation,  and  are  very  accurate.  Their  use 
enables  new  eigenfunctions  to  be  found  by  a  simple  forward  integration  of  the 
differential  equation,  with  little,  if  any,  iteration  necessary. 

Finally,  we  have  studied  the  insulated  wall  problem  using  the  eigen- 

9 

function  series,  a  quadrature  formula  developed  by  Durgin  based  on  Lighthill' s 

4  13 

'  approximation,  an  integral  method  of  Libby  and  Morduchow,  and  an  exact 
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numerical  solution  of  the  partial  differential  equation  obtained  by  the  present 
author.  ^  Comparison  is  made  near  the  initial  plane,  where  the  numerical 
solution  has  a  very  steep  slope  as  it  approaches  the  constant  wall  tempera¬ 
ture  value.  The  Lighthill-Durgin  solution  agrees  almost  perfectly  with  the 
numerical  solution  in  this  region,  having  a  slope  proportional  to  (£/£.- 1)~  '  . 
The'  Fox-Libby  series  differs  considerably  near  the  initial  plane,  even  with 
20  terms,  attesting  to  its  slow  convergence  there.  The  integral  solution  also 
is  inaccurate  there. 

As  £  00  far  downstream,  the  Lighthill-Durgin  solution  and  the  inte¬ 

gral  solution  both  have  the  same  decay  as  the  (presumably  exact)  series 
solution,  but  the  former  is  6%  lower  and  the  latter  8%  higher  than  the  series. 

We  conclude  that  the  Lighthill-Durgin  quadrature  formula  is  an  excel¬ 
lent  approximation  to  the  exact  solution  over  the  whole  range.  This  is  true 
not  only  at  Prandtl  number  unity,  but  also  at  other  values,  since  this  solu¬ 
tion  is  independent  of  that  parameter  when  referred  to  the  adiabatic  wall 

9  10 

temperature.  This  result  is  in  accord  with  Durgin' s  experimental  work,  ’ 
which  was  performed  in  air,  with  Prandtl  number  0.  72.  He  found  the  Lighthill- 
Durgin  solution  to  be  in  very  good  agreement  with  his  laminar  data. 
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IJ  ABSTRACT 

We  study  three  eigenvalue  problems  of  boundary  layer  theory  first 
introduced  by  Fox  and  Libby,  *  and  pursued  further  by  Kotorynski.  ^  By 
completing  Kotorynski's  approximate  analysis,  we  exhibit  approximate 
eigenvalues  and  norms.  With  these,  the  convergence  of  the  eigenfunction 
series  solutions  of  two  heat  transfer  problems  can  be  investigated:  (a)  Flat 
plate  with  one  constant  wall  temperature  region  followed  by  another. 

(b)  Flat  plate  with  a  constant  wall  temperature  region  followed  by  an  in¬ 
sulated  wall.  The  series  for  problem  (b)  is  found  to  converge  everywhere, 
including  at  the  station  of  change  in  the  boundary  condition.  The  series 
for  problem  (a)  does  not  converge  at  that  station,  and  convergence  for  both 
appears  slow  near  that  station.  Then  we  improve  Kotorynski's  approxima¬ 
tion  and  use  the  WKB  method  to  obtain  simple  formulas  for  the  eigenvalues, 
which  prove  very  accurate,  and  enable  the  eigenfunctions  to  be  found  by  a 
simple  forward  integration.  Finally,  we  study  problem  (b)  in  more  detail, 
and  obtain  an  exact  solution  by  a  numerical  finite  difference  procedure. 

This  shows  the  eigenfunction  series  solution  to  be  very  slowly  convergent 
near  the  beginning  of  the  insulated  region.  The  solution  worked  out  by 
Durgin*^  10  based  on  Lighthill's^  approximate  method,  is  very  accurate 
there.  In  addition,  it  is  within  10%  of  the  numerical  solution  everywhere. 

It  is  undoubtedly  the  most  convenient  representation  of  the  solution,  agree¬ 
ing  both  with  the  new  numerical  solution,  and  with  Durgin's  experiments. 
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